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Abstract 

The critical endpoint of the QCD phase diagram is usually expected to belong to the chiral 
critical surface, i.e. the surface of second order transitions bounding the region of first order chiral 
phase transitions for small quark masses in the {m u j,m s ,fi} parameter space. For ft = 0, QCD 
with physical quark masses is known to be an analytic crossover, requiring the region of chiral 
transitions to expand with pt for a critical endpoint to exist. Instead, on coarse N, - 4 lattices, we 
find the area of chiral transitions to shrink with fi, which excludes a chiral critical point for QCD 
at moderate chemical potentials hb < 500 MeV. First results on finer N, = 6 lattices indicate a 
curvature of the critical surface consistent with zero and unchanged conclusions. 



1. Introduction 

The QCD phase diagram has been the subject of intense research over the last ten years. 
Based on asymptotic freedom, one expects at least three different forms of nuclear matter: 
hadronic (low jjg, T), quark gluon plasma (high T) and colour-superconducting (high low 
T), Whether and where these regions are separated by true phase transitions has to be deter- 
mined by first principle calculations and experiments. Since QCD is strongly coupled on scales 
of nuclear matter, Monte Carlo simulations of lattice QCD are presently the only viable approach. 

Unfortunately, the so-called sign problem prohibits straightforward simulations at finite bary- 
on density. There are several ways to circumvent this problem in an approximate way, all of 
them valid for [i/T <1 only lUul]. Within this range, all give quantitatively agreeing results for, 
e.g., the calculation of T c (jj) [3]. Because of the intricate and costly finite size scaling analyses 
involved, determining the order of the transition, and hence the existence of a chiral critical 
point, is a much harder task. Here we discuss the order of the finite temperature phase transition 
as obtained from lattice QCD simulations in the extended parameter space {m u j, m s , T,/i}. 

2. The chiral critical line at ft = 

The schematic situation is depicted in Fig.[TJ beginning with fi = (left). In the limits of zero 
and infinite quark masses (lower left and upper right corners), order parameters corresponding 
to the breaking of a global symmetry can be defined, and one numerically finds first order phase 
transitions at small and large quark masses at some finite temperatures T c (m). On the other hand, 
one observes an analytic crossover at intermediate quark masses, with second order boundary 
lines separating these regions. Both lines have been shown to belong to the Z(2) universality 
class of the 3d Ising model Since the line on the lower left marks the boundary of the 

quark mass region featuring a chiral phase transition, it is referred to as chiral critical line. It has 
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Figure 1 : Left: Schematic phase transition behaviour of Nf = 2 + I QCD for different choices of quark masses at /i = 0. 
On finer lattices, the chiral critical line moves towards smaller quark masses. Right: Critical surface swept by the chiral 
critical line as /i is turned on. Depending on its curvature, a QCD chiral critical point is present or absent, cf. Fig. [2 



recently been mapped out onN, - 4 lattices [ 7]. A convenient observable is the Binder cumulant 
B A (X) = ((X-(X}) 4 }/((X-(X)) 2 ) 2 ,withX = tJ/i{f. At the second order transition, B 4 takes the value 
1.604 dictated by the 3d Ising universality class. In agreement with expectations, the critical line 
steepens when approaching the chiral limit. Assuming the Nf - 2 chiral transition to be in 
the 0(4) universality class implies a tricritical point on the ra s -axis, Fig. Q] (left). The data are 
consistent with tricritical scaling [8] of the critical line with m u d and we estimate m'" c ~ 2.8T C . 
However, this value is extremely cut-off sensitive and likely smaller in the continuum, cf. SecjU 



3. The chiral critical surface 

When a chemical potential is switched on, the chiral critical line sweeps out a surface, as 
shown in Fig. [T] (right). According to standard expectations [8], for small m u j, the critical line 
should continuously shift with // to larger quark masses until it passes through the physical point 
at /i£, corresponding to the endpoint of the QCD phase diagram. This is depicted in Fig.Q](right), 
where the critical point is part of the chiral critical surface. However, it is also possible for the 
chiral critical surface to bend towards smaller quark masses, cf. Fig.|2](right), in which case there 
would be no chiral critical point or phase transition at moderate densities. For definiteness, let 
us consider three degenerate quarks, represented by the diagonal in the quark mass plane. The 
critical quark mass corresponding to the boundary point has an expansion 



— — = 1 + > Q — . (1) 
m c (0) f-f KnTJ 

A strategy to learn about the chiral critical surface is now to tune the quark mass to m c (0) and 
evaluate the leading coefficients of this expansion. In particular, the sign of c\ will tell us which 
of the scenarios is realised. 

The curvature of the critical surface in lattice units is directly related to the behaviour of the 
Binder cumulant via the chain rule, 
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Figure 2: /^-dependence of the Binder cumulant on the chiral critical line for Nf = 3 (left) [ 9||. Thus, the scenario on the 
right is realised. For heavy quarks the curvature has also been determined | 6] and the first order region shrinks with fi. 



While the second factor is sizeable and easy to evaluate, the //-dependence of the cumulant is 
excessively weak and requires enormous statistics to extract. In order to guard against systematic 
errors, this derivative has been evaluated in two independent ways. One is to fit the corresponding 
Taylor series of B4 in powers of [i/T to data generated at imaginary chemical potential Blah the 
other to compute the derivative directly and without fitting via the finite difference quotient (0] 

8B 4 B 4 (afi)-B 4 (0) 

= hm — . (3) 
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Because the required shift in the couplings is very small, it is adequate and safe to use the original 
Monte Carlo ensemble for am c Q ,[i = and reweight the results by the standard Ferrenberg- 
Swendsen method. Moreover, by reweighting to imaginary fi the reweighting factors remain real 
positive and close to 1 . The results of these two procedures based on 20 and 5 million trajectories 
on 8 3 x4, respectively, is shown in Fig.[2](left). The error band represents the first coefficient from 
fits to imaginary fi data, while the data points represent the finite difference quotient extrapolated 
to zero. Both results are consistent, and the slope permits and extraction of the subleading // 
coefficient. After continuum conversion the result for Nf — 3 is c\ = -3.3(3), C2 = -47(20) |9;]. 
The same behaviour is found for non-degenerate quark masses. Tuning the strange quark mass 
to its physical value, we calculated ni c ' d (p) with c\ = -39(8) and C2 < 0, Fig. [3] (left). Hence, on 
coarse N, - A lattices, the region of chiral phase transitions shrinks as a real chemical potential is 
turned on, and there is no chiral critical point for /j^ <500 MeV, as in Fig. [2] (right). Note that one 
also observes a weakening of the phase transition with // in the heavy quark case |g], in recent 



model studies of the light quark regime 11121 11311 . as well as a weakening of the transition with 
isospin chemical potential B14ll . 



4. Towards the continuum, N t = 6 

The largest uncertainty in these calculations by far is due to the coarse lattice spacing a ~ 0.3 
fm onN, = A lattices. First steps towards the continuum are currently being taken on N, = 6, a ~ 
0.2 fm. At /j. = 0, the chiral critical line is found to recede strongly with decreasing lattice spacing 
fnifrill : for Nf = 3, on the critical point m„{N t = A)lm Jt {N t — 6) ~ 1.8. Thus, in the continuum 
the gap between the physical point and the chiral critical line is much wider than on coarse 
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Figure 3: Left: ^-dependence of the Binder cumulant on the chiral critical line for Nf = 2 + 1 on N, = 4 with physical 
strange quark mass. Right: The same for Nf = 3, N, = 6. 



lattices, as indicated in Fig. Q] (left). Preliminary results for the curvature of the critical surface, 
Fig.[3](right), result in c\ = 7(14), -17(18) for a LO,NLO extrapolation in /j 2 , respectively. Thus 
the sign of the curvature is not yet constrained. But even if positive, its absolute size is too small 
to make up for the shift of the chiral critical line towards smaller quark masses, and one would 
again conclude for no chiral critical point below fig <500 MeV in this approximation. Higher 
order terms with large coefficients would be needed to change this picture. 

However, on current lattices cut-off effects appear to be larger than finite density effects, 
hence definite conclusions for continuum physics cannot yet be drawn. A general finding is 
the steepness of the critical surface, making a possible critical endpoint etremely quark mass 
sensitive. Furthermore, the shrinking of the critical quark masses with diminishing lattice spacing 
makes it likely that the tricritical point moves to m u j = 0, m f " e < m p s ys , in which case a possible 
critical endpoint might belong to an entirely different critical surface. 



Acknowledgments 

This work is partially supported by the German BMBF, project No. 06MS254. 



References 

[1] O. Philipsen, Eur. Ph ys. J. ST 152 (2007) 29 larXiv:0708.1293l [hep-lat]]. PoS LAT2005 (2006) 016 

jarXiv:hep-lat/0510077]. 

[2] C. Schmidt, PoS LAT2006 (2006) 021 |arXiv:hep-lat/0610116|. 

[3] P. de Forcrand and S. Kratochvila, PoS LAT200S (2006) 167 |hep-lat/0509143 |. 

[4] F. Karsch, E. Laermann and C. Schmidt, Phys. Lett. B 520 (2001) 41 | arXiv:hep-lat/0107020 |. 

[5] P. de Forcrand and O. Philipsen, Nucl. Phys. B 673 (2003) 170 | arXiv:hep-lat/0307020 |. 

[6] S. Kim, Ph. de Forcrand, S. Kratochvila and T. Takaishi, PoS LAT2005, (2006) 166 | arXiv:hep-lat/05 10069 |. 

[7] P. de Forcrand and O. Philipsen, JHEP 0701 (2007) 077 [hep-lat/0607017 |. 

[8] A. M. Halasz et aL, Phys. Rev. D 58 (1998) 096007 | arXiv:hep-ph/9804290 1. 

[9] P. de Forcrand and O. Philipsen, JHEP 0811 (2008) 012 | arXiv:0808.1096 [hep-lat]]. 

[10] P. de Forcrand, S. Kim and O. Philipsen, PoS LAT2007 (2007) 178 | arXiv:07 11.0262 [hep-lat]]. 

[11] G. Endrodief a/.,PoS LAT2007 (2007) 182 |arXiv:07 10.0998 [hep-lat]]. 

[12] K. Fukushima, Phys. Rev. D 78 (2008) 1 14019 | arXiv:0809.3080 [hep-ph]L 

[13] E. S. Bowman and J. I. Kapusta, Phys. Rev. C 79 (2009) 015202 | arXiv:0810.0042 [nucl-th]]. 

[14] J. B. Kogut and D. K. Sinclair, Phys. Rev. D 77 (2008) 114503 | arXiv:0712.2625 [hep-lat]]. 



4 



